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Answers to some salient questions, which arise in quantum plasmas, are given. Starting from the 
Schrodinger equation for a single particle it is demonstrated how the Wigner-Moyal equation can 
be derived. It is shown that the Wigner-Moyal type of equation also exists in the classical field 
theory. As an example, from the Maxwell equations the Wigner-Moyal type of equation is obtained 
for a dense photon gas, which is classical, concluding that the Wigner-Moyal type of equation can 
be derived for any system, classical or quantum. A new type of quantum kinetic equations are 
presented. These novel kinetic equations allow to obtain a set of quantum hydrodynamic equations, 
which is impossible to derive by the Wigner-Moyal equation. The propagation of small perturbations 
^ ' and instabilities of these perturbations are then discussed, presenting new modes of quantum plasma 

waves. In the case of low frequency oscillations with ions, a new Bogolyubov type of spectrum is 
found. Furthermore, the Korteweg-de Vries (KdV) equation is derived and the contribution of the 
Madelung term in the formation of the KdV solitons is discussed. 
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I. INTRODUCTION 

Our understanding of the collective behavior of quantum plasmas has recently undergone some appreciable theoret- 
ical progress [1] , [2] , [3] . Quantum plasma is of primary interest in its potential application in modern technology, such 
as metallic and semiconductor nanostructures. As well as, the quantum plasmas are prevalent in planetary interiors 
Oh, and in compact astrophysical objects, e.g., the interior of white dwarf stars, magnetospheres of neutron stars and 
ryj ■ magnetars, etc. Despite extensive theoretical efforts in the remote past [4]- [8], where the properties of linear electron 
oscillations in a dense Fermi plasma have been studied, and more recent ones [9], [10], there are questions and issues 
c/j which has to be clarified. 

>V The purpose of the present work is to elucidate and give the correct formulation to quantum plasmas, and demon- 
strate which equations can be used, as well as which are invalid for the description of the collective behavior of 

i quantum plasmas. So that in this paper, we provide answers to some salient questions, which arise in quantum 

plasmas. We start from the Schrodinger equation for a single particle and demonstrate how the Wigner-Moyal equa- 

t-H ■ tion can be derived. We then show that the Wigner-Moyal type of equation also exists in the classical field theory. 

Next we present a new type of quantum kinetic equations. These novel kinetic equations allow us to obtain a set of 
quantum hydrodynamic equations, which is impossible to derive by the Wigner-Moyal equation. We then discuss the 
propagation of small perturbations and instabilities of these perturbations, presenting new modes of quantum plasma 
waves. Finally, we derive the Korteweg-de Vries (KdV) equation in order to discuss the solitary waves in quantum 
♦ plasmas. 



, — i 
>< 



II. DERIVATION OF THE WIGNER-MOYAL EQUATION 

We start with a single particle and for this purpose, we employ the non-relativistic Schrodinger equation 



B^(r t) ti 2 

<3 ; lh ^T^ + 2^ A * (f ' * } " v(f > w> t] = ' (1) 

where t) is not only the wave function of the single particle, but also an amplitude of the probability. Namely, 
| ^(r, t) | 2 df is the probability of finding values of coordinates of the particle within the interval f, r + dr. The 
probability density is | >F(r, t) | 2 , and | <&(f, t) \ 2 ~ k. 

Note that the Schrodinger equation (1) describes the wave properties of the single particle. For the dense system, 
N-body (N — > oo) system with an interaction potential V(r N ) = V{r\,?2, ...fjv) the Schrodinger equation reads 



dV(r N ,t) 



h 2 N d 2 

k=l k 



^{r l \t). (2) 



Obviously, in this case it is very difficult to discuss the wave properties. Therefore, we should proceed from waves to 
particles. 
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We consider Eq.(l) at two different points, r\ and r 2 

in ' + ^Ax^(fi,t) - V{r u t)Wut) = , (3) 

d^*(fo t) h 2 

- ih—^ 1 + ^A 2 **(r 2 ,i) - nr 2 ,t)**(f 2l t) = . (4) 
We now multiply Eq.(3) by ^*(r 2 ,i) and Eq.(4) by \&(ri,i), and subtract the resulting equations to obtain 

*(fi,t)**(f?,,t) + _- (Ai- A 2 )^(r 1 ,t)^*(f2,t)-[V(r 1 ,t)-V(f2,t)}^(r 1 ,t)^%f 2 ,t)=0 . (5) 



Introducing a new variables p = f\ — f 2 , r = ^(fi + r 2 ), so that 

d d \ ( d d 



A A — ( ]( 1—2 

1 2 V 9r*i <9r 2 ' \ dr\ dr 2 ' 



drdp 



*(n,t)**(f 2 ,t) = *(r + |,i)**(r- |, t) , 



and multiplying Eq.(5) by e , and integrating over p we get 



(6) 



where 



/ = U = j dpe-™ W+ P -,t) (7) 

is the Wigner function. 

In order to derive the Wigner-Moyal equation one should suppose that V(r, t) is a smooth function. Expansion 
v(r + f, tj — v(f — %,tj in Taylor series leads to the Wigner-Moyal equation 

^+ {1 ^f™-l sin & F -V?) V(f,t)f w =0. (8) 

It should be emphasized that this equation has the following restrictions: V(f, t) should be the smooth function, and 
Eq.(8) is not suitable for the construction of fluid equations, even if we apply it to a macroscopical body. 
Note that the integration of the Wigner function over momenta gives the configurational probability density 
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f w dp=\*(r,t)\ 2 , (9) 

whereas the integration over spatial coordinates gives the momentum probability 

f w df=\ <p[p,t) | 2 . (10) 

Moreover, the mean value of any physical quantity L(p, q) can be calculated by the statistical distribution function 

<L>= J L{p,q)f w (f,p)dfdp . (11) 

However, the Feynman notes [11] "Although Wigner function f w (f,p) satisfies Eqs.(9) and (10), it cannot be regarded 
as the probability for finding the particle at the point r and with the momentum p, because f w (r , p) can become 
negative for some values of p and r. Also, Eq.(ll) is not true for a general function L(p,q)". Thus the Wigner 
distribution function has no physical meaning. 
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We here also write the Wigner distribution function for N-particle system 

The question is: does the Wigner's type of equation exist in the classical approximation, or it is purely quantum? 
We emphasize that such type of equation one can derive for any system, classical or quantum. To show this, we can 
employ the relativistic Maxwells equation 

. t* d 2 A n - , r 

ot 7 V m^c 4 

and consider it at two distinct points (f i,r2). Following the procedure described above for the Wigner-Moyal equation, 
one can obtain for the occupation number N(f,k,t) of photons [12], [13], [14] 

J^iV(r, k, t) + (v g • V)7V(r, k, t) - Lo 2 p sin(^V P • V E ) V{r, t)N(f, k,t)=0 (13) 

and 

/dk 
j^N{r,k,t) » 1 , 

which describes the dense photon gas. 



III. NEW KINETIC EQUATIONS 



Some uncertainty of the Wigner-Moyal equation demands the derivation of a new quantum kinetic equation incor- 
porating all the quantum properties of the body. Moreover, from which must follow a set of fluid equations with a 
quantum term. To this end, we start with a single fermi particle and employ the non-relativistic Pauli equation [15], 
which reads 



in—^-. — h 



dt 



2m„ 



-A*„ - 



ieh . t „ _ ? , e 2 A 2 
-(A-V + VA) ■ 



2m n c 



2m n c 2 



e a cp - fl a ■ H 



, 



(14) 



where $f a — $> a (r,t,ij) is the wave function of the single particle species a, having the spin s = 1/2(7 (a = ±1). 
A(r,t) and f(r,t) are the vector and scalar potentials, respectively. The last term in Eq.(14) is the potential energy 
of magnetic dipolc in the external magnetic field, the magnetic moment of which is 



Ma = 



eh 
2m„c 



a = n p a 



(15) 



where fj,p is the Bohr magneton and a is the operator of the single particle [15], [16]. 
Use of the Madelung representation [17] of the complex function \I> Q 



* a (f, t, a) = a a (f, t, a) exp 



iS a (f,t,a) 



where a a (f,t,a) and S a (f,t,a) are real, in the Pauli equation (14), yields the following two equations 

dal(f,t,a) 



dt 



+ V.(otat,*)£)=0 



(16) 



(17) 



d ^ = e (E+ ^JL) + ilviv 2 a Q + W V(a • H) . (18) 

Equation (17) has an obvious physical meaning. Namely, a 2 a = \^ a \ 2 is the probability density of finding the single 
particle at some point in space with a spin s. Whereas, p a — VS* — ^rA is the momentum operator of the particle. 
Note that if there is no spin dependence of the wave functions ^> a (f, t), then p a (r, t) becomes the ordinary momentum 
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of the particle. It should be emphasized that in Eq.(18) the second, quantum Madelung, term describes the diffraction 
pattern of a single electron [18]. 

In the case when there are no external electric and magnetic fields (E = 0, H = 0), after the linearization of Eqs. 
(17) and (18), we get the frequency of quantum oscillations of a free electron 

Based on the diffraction pattern of electrons, Born has given a statistical interpretation of the wave function, 
which states that in every point in space at a given time the intensity of the de Broglie waves is proportional to the 
probability of observing a particle at that point in space. Therefore, a density of probability distribution fs in phase 
space for the single particle is 

l*| 2 = J d 3 P f s (f,p,t) . (20) 

This function fs, obviously must satisfy the normalization condition over all phase space 

J d 3 r J d 3 p fs = J d 3 r |*| 2 = 1 . (21) 

Moreover, Bogolyubov has introduced a one particle distribution function [19] for the system as a whole from Liouville's 
theorem regarding the distribution function f„(t, t\, t%...tn), (where T a is the set of coordinates and momentum 
components for the a particle), and derived the Vlasov and Boltzmann equations in the gas approximation, which 
means that the plasma parameter rj (representing the ratio of the average potential energy < U > of particles 
interaction to the average kinetic energy < Ek >) must be less than unity, i.e., r\ = ^3 <C 1. Note that the one 
particle distribution function f{r,p,t) is normalized to unity, whereas the Liouville's function f N to total number of 
particles, i.e., f N — Nf(f,p,t). The same relation between f N and f(f,p,t) in an alternative description of kinetic 
theory was obtained by Klimontovitch [20], [21]. To make it more lucid, we shall give a simple explanation about a 
single particle and one particle distribution function [21], [16]. Namely, the probability density of the single particle is 
one particle per unit volume, n s (f, t) = \^(f, t)\ 2 , with dimensions 1/V. Whereas the one particle distribution function 
means that in spite of the large number of particles in the unit volume all of them have only one r and p. This permits 
us to express the number density of particles per unit volume as n(f,t) — J d 3 p f(f,p,t) = N/V. Therefore, we can 
write 

n(r, t) = Nn s (r, t) = N\*(F, t)\ 2 = J d 3 p f(r,p, t) . (22) 

Thus, n(r, t) is the density of quantum particles per unit volume. We have assumed that the total number of particles 
of each species is conserved. 

Non-equilibrium states of a Fermi quantum gas are described by the one particle distribution function f a (r , p, t, a) , 
which satisfies the quantum Boltzmann equation 

df a (f,p,t,a) , ^, 7 \ ^ s) + d& W, t, gj = C{L) , (23) 



dt \m ) ' ' dt dp 

Equation (23) for quasi-particles in a Fermi liquid was written by Landau [22], [23]. Here C(f a ) is the collision integral, 
which describes the variation of the distribution function due to particle collisions, the derivative is determined 
by the force acting on the particle, the expression of which is given by the equation (18). When the spin of particles 
is taken into account, the distribution function f a becomes an operator with respect to the spin variables a. In this 
case the total number density of particles n a (r, t) equals 

n a (r, *) = E / f° ^ & *' *) • ( 24 ) 

If there is no spin dependence of the distribution function, then f a becomes the ordinary quasi-classical distribution 
function f a (r,p,t). Note that the condition for quasi-classical motion is that the de Broglie wavelength Xb — fi/pF 
of the particle must be very small compared with the characteristic length L, over which f(r,p, t) varies considerably. 

If collisions between particles were entirely negligible, each particle of the system would constitute a closed sub- 
system, i.e., one could neglect the collision integral in Eq.(23), and the distribution function of particles would obey 
Liouville- Vlasov equation. 
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Hereafter, we consider the system as a spinless. We substitute the equation of motion of single particle (18) 
(neglecting the last term) into the kinetic equation (23) taking into account the definition of the density of particles 
(22) to obtain 

— + (v-V)f a +e a [E + — + V^=Ay^ — = C(f a ) . (25) 

at V c / op 2m a ^fn^ Op 

It should be emphasized that this is a novel equation with the quantum term, which contains all the information on 
the quantum effects. We specifically note also that this equation is rather simple from the mathematically point of 
view. 



IV. PROPER WAVES IN A FERMI QUANTUM PLASMA 

We now employ the equation (25) to study the propagation of small longitudinal perturbations (H = 0, E = —Vip) 
in an electron-ion collisionless plasmas. For a weak field, we look for the electron and ion distribution functions in 
the form f a — f a o + 5f a , where f a o is the stationary isotropic homogeneous distribution function unperturbed by the 
field, and 5f a is the small variation in it due to field. After linearization of Eq.(25) with respect to the perturbation, 
we assume Sf a and Sip vary like expi(k ■ r — uit). 



Using the Poisson's equation 



and assuming the Fermi degeneracy temperature Tp = jf^ (Kb is the Boltzmann coefficient, the Fermi distribution 

2 

function is the step function f a o — Q(e F a — e), where e Fa = m °2 F ° ) much more than the Fermi gas temperature, 
then we obtain after some algebra the quantum dispersion equation 

3Wp Q f UJ UJ + kVFa \ 



where 

r -l i 3 ^ fc2 (i - \ n UJ + kvFa ) 

Am 2 a v F a V 2kv Fa u-kv Fa > 

and u> can be more or less than kvFe- Note that for to ^> kvpi, ^ w 1. 

Let us first consider the electron Langmuir waves, supposing that the ion mass m, — > oo and w » kvp e , or the 
range of fast waves, when the phase velocity exceeds the Fermi velocity of electrons. In this case we get the dispersion 
relation 

- 2 =< + ^ + ^ + - (28) 

which has been previously derived by Klimontovich and Silin [5]. The expression (28) exhibits that the high frequency 
oscillations of electrons of a degenerate plasma remain undamped in the absence of particles collisions. Note that the 
Landau damping is also absent, since according to the Fermi distribution there are no particles with velocities greater 
than the Fermi velocity which could contribute to the absorption. 

As k increases, Eq.(28) becomes invalid, but still uj > kvp e and the Landau damping is absent. We now introduce 

the Thomas- Fermi screening wave vector krF = In the limit fc 2 k\ F , w tends to kvpe (at m, — > oo) and 

we obtain from (27) 

2(_|L + J^L) 

uj = kv Fe (l + 2 cxp{ kTF J" }) . (29) 

pe 

If we neglect the quantum term uj q in Eq.(29), then we recover waves known as the zero sound, which are the 
continuation of the electron Langmuir wave (28) into the range of short wavelength. Thus the expression (29) 
represents the quantum correction to the zero sound. 



6 



Special and very important case in a quantum plasma is an one-fluid approximation. In this case one assumes 
that the characteristic dimension R of inhomogencitics in the plasma is larger than the electron Thomas-Fermi length 
\ TF = . In the perturbed plasma the potential of electric field is determined by the Poisson's equation 

V OUJp e 

Aip = 47re(<5n e — ZiSrii) . (30) 

To estimate a magnitude of the term on the right-hand side of Eq.(30), we consider two cases e<p <C £f and ap <~ Ef- 
Noting that Atp <~ f/R 2 , for the weak electric field Sn e = f^ofj we rewrite Eq.(30) in the form 

I fae - ZjSrij \ 2 TF 

I ^ I- "jjT « 1 ' V L > 

This inequality (31) is also satisfied for the strong perturbation, eip ~ Ef and 5n e ~ n e . Thus we can conclude 
that the uncompensated charge density is small compared to the perturbation of the electron and ion charge density 
separately. So, for Xtf -C R we further assume that the quasi neutrality 

n e = m (32) 

is satisfied. This equation along with the equation of motion of ions and the equation giving the adiabatic distribution 
of electrons allow us to define the potential field. In such approximation the charge is completely eliminated from the 
equations, and the Thomas-Fermi length Xtf disappears with it. 

In order to construct the one-fluid quantum kinetic equation, we neglect the time derivative in the equation (25) 
of electrons, as well as the collision terms, suppose E — —Vip and H = 0, and write the dynamic equations for the 
quasi-neutral plasma (32) 

(„. V )/. + v(e», + |--LA^)f =0, (33) 



^ + („.V)/,-Ve v ^=0. (34) 
In Eq. (34) we have neglected the quantum term as a small one. Note that the Fermi distribution function of electrons 

fe = 7~^2 - S (35) 

exp J ^-M + l 



satisfies the equation (33). Here U = etp+ -^-^Ay/n , and jj, e is the chemical potential. 

For the strongly degenerate electrons, i.e., T e — > (fi e = Ef), the Fermi distribution function becomes the step 
function 

f e = 6(e F + U- £-) , (36) 
which allows us to define the density of electrons (n e = n, = n = 2 j r^jsf) 
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We now express eip from the equation (37) and substitute it into the kinetic equation (34) to obtain 



This is the nonlinear kinetic equation of quantum plasma in the one-fluid approximation, which incorporates the 
potential energy due to degeneracy of the plasma and the Madelung potential. 
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To consider the propagation of small perturbations / = fo(p) + 8f{r,p,t) and n = uq + Sn(f,t), we shall linearize 
Eq.(38) with respect to the perturbations, look for a plane wave solution as e l ( fc r ~"*\ and derive the dispersion 
equation, which resembles the Bogolyubov's dispersion relation in the frequency range kvFe ^> w » kvp i} 



P 2 F 



Ti 2 k 2 , 2e Fe h 2 k 2 

= k \hz- + 7zrz^- (39) 



y 3m e rrii 4mim e y Zrrii 4m,im e 

We specifically note here that this type of spectrum (39) was derived by Bogolyubov for the elementary excitations 
in a quantum Bose liquid, and the microscopic theory of the superfluidity of liquid helium was created [23], [24]. 

How one can explain the similar spectrum (39) in the Fermi gas? In the dense and low temperature plasma the 
formation of bound states is possible due to the attractive character of the Coulomb force [16], [25], [26], [27]. As is well 
known, when the density of particles increases and the temperature goes to zero, the nuclear reaction leads to the 
capture of electrons by nuclei. In such reaction the charge on the ions (nucleus) decreases. Because we have assumed 
the quasi-neutrality, we have therefore supposed that all electrons are in the bound state with ions. Thus one can say 
that in the one-fluid approximation the Fermi plasma may become the Bose system due to the bound state or this 
approximation may imply the formation of the Bose atoms. 

The same Bogolyubov's type of spectrum (39) follows from Eq.(27) in the range of intermediate phase velocities 
kvFe S> cj >• kvpi- In this case, except the expression for the real part of frequency (39), we get the imaginary part 
of ijj from the dispersion relation (27) 

ImL0 = -^k^, (40) 
12 m, v ' 

which is much less than Reoj. The expression (40) indicates that the electrons play role in the absorption of oscillations, 
since their random velocities greatly exceed the phase velocity. So that the damping rate is determined by the electrons 
alone. 

V. DERIVATION OF QUANTUM HYDRODYNAMIC EQUATIONS 

We shall now derive a set of fluid equations. The Boltzmann and Vlasov type of quantum kinetic equations (25) 
and (38) give a microscopic description of the way in which the state of the plasma varies with time. It is also well 
known how the kinetic equation can be converted into the usual equations of fluids. Following the standard method, 
we can derive the equations of continuity and motion of macroscopic quantities from Eq.(38) 

rim 

™ + V • (nit) = (41) 

— + (t»-V)tr= 2_^v(— 2/3 + ^ V — A^, 42 

ot rrii n 2m e m,i sjn 



where u(r,t) is the macroscopic velocity of the plasma 



Obviously from Eqs.(41) and (42) after linearization follows the Bogolyubov's type of dispersion equation (39). 

The more general set of fluid equations for a kind of particles we can obtain from the equation (25) taking into 
account that any elastic scattering should fulfill the general conservation laws of particle number, momentum and 
energy 

^+V-(«A)=0, (44) 
d<p a > _ _ ^. 1^ fA 1 „_ h 2 „ 1 



+ (u a • V) <p a >= e a (e + -u a x B] — VP Q + V— — A x /n^ 

V c / n„ 2m„ y/n^ 

Pa C(f a ) , (45) 



dt 

1 f 2dp 



n a J (2irh) 3 
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where P a - ^ / 



^}+i • 

Moreover, the equation of the internal energy of each particle species can be derived from the kinetic equation (25) 

2 

by multiplying it by , integrating over the momentum and employing the equations of continuity (44) and motion 
(45). In the non-relativistic limit the result is 

d i < £a > 1 f 2d3 P ni r \ (ar\ 

dt x *^nr = ^rj-; J e « c ^ ' (46) 

where e a — ^ Pa ~^ a> ^ and < e a >= J (f^p £ a fa is the internal energy of particles, and the collision terms can 
be positive or negative. 



VI. DERIVATION OF THE KORTEWEG-DE VRIES EQUATION 

We now consider the propagation of weakly nonlinear ion-sound waves. To this end, we generalize the Bogolyubov's 
type of dispersion relation (39), where we assumed the quasi- neutrality (32), taking into account the small separation 
of charges. For the linear ion-sound waves (kvp e »u» kvpi) we have 

lo = kv s + /?fc 3 - ak 3 , (47) 



where v s = \J m e /3m,iVF e is the ion sound velocity, f3 = |A^z; s , and a — v s X?n F /3 is due to the separation of 
charges. 

Equation (47) describes sufficiently long waves in plasmas, where cj/k at k — > has a finite limit (weakly dispersive 
waves) . The nonlinear waves for the weakly dispersive plasmas are governed by the Korteweg-de Vries equation, which 
we shall derive in the following. 

As is well known the nonlinearity, generally speaking, causes the distortion of wave profile, which steepens until the 
dissipative effects come into play. The later causes the broadening of the profile and finally balances the nonlinear 
steepening, so that the nonlinear stationary waves can be formed in the dispersive plasma. To show this, we follow 
the standard method [28]- [30], and from Eq.(37) obtain the expression for the density of electrons 

3 /n\ 2 / 3 

n e = n i + -n Q \ 2 TF A(^) , (48) 
2 vno / 

where n = n 0e = nm is the unperturbed density. Note that the second term in Eq.(48) is small since A^ F <§C R 2 , 
where R is the characteristic scale length of the perturbation. 
In one dimensional case equations (41) and (42) now read 

<9n, d 



du d EFe d /n e \ 2 / 3 h 2 d 1 

dt dx mi dx\no) 2m e mi dx yJrT e 

Substitution of Eq.(48) into Eq.(50) yields the equation of motion for ions 

du . d 3 2 d (, . d 2 \/n,\ 2 / 3 . h 2 did 2 



RT e . (50) 



dt ^ dx 2 s dx \ TF dx 2 ) \no) 2m e mi dx y/ni dx 2 ^ ^ 

Here the last two terms describe the dispersive effects. 

Obviously from Eqs.(49) and (51), which formally resemble the Boussinesq equations, after the linearization follows 
the dispersion equation (47). If we neglect the dispersive terms in Eq.(51), then Eqs.(49) and (51) reduce to the Eulcr 
equations. 

We now consider the nonlinear one dimensional traveling waves in a single direction, so-called simple waves. Since 
the velocity u and the density rn in such a wave are functions of the same variable x — ut, we suppose that for a wave 
with any amplitude, all quantities describing a simple wave are represented in the form of functions of one of these 
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quantities, i.e. it can be assumed that u — u(n) or n = n(u). From Eqs.(51) (without the dispersive terms) and (49) 
then follows for a wave propagating in the positive x direction 



du 
drii 



1/3 2/3 
n n i 



or 



and 



t 



—u 

3 



3v, 



F(u) 



(-) 

Vn / 



1/3 



1 



(52) 



(53) 



where F(u) is an arbitrary function of the velocity. 

Expressions (52) and (53) are general solutions of arbitrary traveling waves and determine the velocity as an implicit 
function of x and t. 

We can also rewrite Eq.(53) in the form 



v(u) =v s + -u , 



(54) 



u = Q(x — v(u)t) , 

where Q(x) is an arbitrary function defining the wave shape at t=0 

Q(x)=u(x,0) . 

Since the velocity u is the function of density, therefore it is different for different points at the wave profile, i.e., the 
profile changes in the course of time. One can see from Eq.(52) that du/dn > 0, i.e., the velocity of propagation of a 
given point at the wave profile increases with the density. This condition du/dn > is of the shock waves [28]- [30]. 
At the point, where the shock is formed, the dispersive terms that we have neglected in Eq.(51) become important 
because the nonlinear steepening can be compensated by the dispersive spread. For some conditions the formation of 
solitary waves takes place. 

Expressing the ion density through the velocity from Eq.(52) and substituting it into Eq.(51), we obtain KdV 
equation 



du 

~dt 



+ 



4 \du d 3 u 

+ r)d- x =- {A - B) ^ 



(55) 



where A = v s \^, F and B — h /4m e miV s . 

We emphasize here that Eq.(55) has two physically distinct solutions. Namely, if the right-hand side of this equation 
is negative, then the solution of Eq.(55) corresponds to the compressional soliton. Whereas, in the opposite case, 
when the quantum (Madelung) term exceeds the term due to the charge separation, the solitary wave is a rarefaction 
wave. 

Interesting case is when A « B. In this case there is no more balance between the terms of the charge separation 
and the quantum, and the nonstationary case of shock formation takes place. The condition A w B is satisfied for 
the density n <~ 10 25 cm~ 3 . 



VII. SUMMARY 



Starting with the Schrodinger equation for a single particle, we have derived the Wigner-Moyal equation and 
pointed out the limitations of the usability of this equation. Namely, we have revealed that the Wigner quantum 
kinetic equation is not appropriate for the construction of fluid equations. It is elucidated that the Wigner distribution 
function has no physical meaning. We have also demonstrated that the Wigner-Moyal type of equation can be derived 
for any system, classical or quantum. We have obtained a new type of quantum kinetic equation motivated by its 
necessity for the correct description of the collective behavior of the Fermi quantum plasma. This novel quantum 
kinetic equation allows to obtain a set of quantum hydrodynamic equations. In our formulation of kinetic, as well as 
hydrodynamic equations the quantum term, which contains all the information on the quantum effects, is incorporated. 
We have study the propagation of small longitudinal perturbations in an electron-ion collisionless plasmas, deriving a 
quantum dispersion equation. We have discussed the quantum correction to the zero sound. For the special interest 
case we have derived a quantum kinetic equation in the one-fluid approximation and have shown that, when a charge 
and the Thomas-Fermi screening length are not present in the equation, the solution of linear kinetic equation leads 
to the Bogolyubov's type of dispersion equation, which is valid in the Bose fluid. It is clear that in the dense strongly 
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coupled plasmas with attractive interaction the formation of the bound state is possible, so that a neutral Bose atoms 
can be created. Furthermore, we have derived KdV equation in order to discuss the nonlinear effects in a quantum 
plasma, and found that due to the quantum term the compressional solitons may become rarefaction waves. These 
investigations may play an essential role for the description of complex phenomena that appear in dense astrophysical 
objects, as well as in the next generation intense laser-solid density plasma experiments. 
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